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Abstract 

The presence of signal outage, due to shadowing and blockage, is expected to be the main bottleneck 
in millimeter wave (mmWave) networks. Moreover, with the anticipated vision that mmWave networks 
would have a dense deployment of base stations, interference from strong line-of-sight base stations 
increases too, thus further increasing the probability of outage. To address the issue of reducing 
outage, this paper explores the possibility of base station cooperation in the downlink of a mmWave 
heterogenous network. The main focus of this work is showing that, in a stochastic geometry framework, 
cooperation from randomly located base stations decreases outage probability. With the presumed vision 
that less severe fading will be experienced due to highly directional transmissions, one might expect 
that cooperation would increase the coverage probability; our numerical examples suggest that is in 
fact the case. Coverage probabilities are derived accounting for: different fading distributions, antenna 
directionality and blockage. Numerical results suggest that coverage with base station cooperation 
in dense mmWave systems and with no small scale fading considerably exceeds coverage with no 
cooperation. In contrast, an insignificant increase is reported when mmWave networks are less dense 
with a high probability of signal blockage and with Rayleigh fading. 
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I. Introduction 

One of the fundamental goals for 5G is a radieal inerease in data rates Q. It is antiei- 
pated that higher data rates will be aehieved by extreme densifieation of base stations, massive 
multiple-input-multiple-output (MIMO), increased data rate and/or base station cooperation [j^. 
However, prime microwave wireless spectrum has become severely limited, with little unassigned 
bandwidth available for emerging wireless products and services. Therefore, to fulfill the need 
for increased bandwidth, millimeter wave (mmWave) spectrum between 30 and 300 GHz have 
been considered for future 5G wireless mobile networks. Until recently, mmWave frequency 
bands were presumed to be unreliable for cellular communication due to blockage, absorption, 
diffraction, and penetration, resulting in outages and unreliable cellular communications 0 . 
However, the advances in CMOS radio-frequency circuits, along with the very small wavelength 
of mmWave signals, allows for the packing of large antenna arrays at both the transmit and 
receive ends, thus providing highly directional beam forming gains and acceptable signal-to-noise 
ratio (SNR) Q, Q. This directionality will also lead to reduced interference when compared to 
microwave networks [[^. It is thus anticipated that mmWave spectrum holds tremendous potential 
for increasing spectral efficiency in upcoming cellular systems 0 . 

To further address the demand for higher data rates, cooperation between macro, pico and 
femto base stations has been proposed to enable a uniform broadband user experience across 
the network. The dynamic coordination across several base stations - known as coordinated 
multipoint (CoMP) - will limit the intercell interference thus increasing throughput and enhancing 
performance at cell borders [|^. 

Past Work: Coverage and capacity in mmWave cellular systems and in CoMP networks have 
been studied. In [|^ the authors compared the performance, in terms of coverage and capacity, of 
a stochastic geometry based mmWave network (without CoMP) to a microwave cellular network, 
at a single antenna receiver (typical user). In [|^, directionality at the transmitters, intra-cell and 
inter-cell interference were accounted for but blockage was not included in the analysis. The 
authors show that coverage in mmWave systems increases with the decrease in the half-power 
beam width of the radiation pattern. In fact, having narrower beams decreases beam overlap, 
thus decreasing intra-cell and inter-cell interference and increasing coverage probability. In this 
paper, we propose to study the problem of base station cooperation in the downlink of dense 
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mmWave heterogenous network as a means to eombat bloekage and deerease signal outage. Our 
derivations of the eoverage probability, similarly to Q, aeeount for interferenee experieneed at 
the typieal user, but in addition bloekage is ineorporated in the analysis. 

In [|7| (see also journal version in Q) the authors proposed a stoehastie geometry framework 
to evaluate the performanee of mm Wave eellular networks (without CoMP) with bloekage. The 
authors ineorporate bloekage by modeling the probability of a eommunieation link - being either 
a line-of-sight (LOS) or non-LOS (NLOS) link - as funetion of the length of the eommunieation 
link from the serving base station. Different pathloss laws were applied to the LOS and NLOS 
links. Numerieal results in 0 suggest that higher data rates ean be aehieved when eompared to 
mierowave eellular networks. One of the interesting observations made in Q is that mmWave 
networks should be dense but not too dense - sinee the number of LOS interfering base stations 
inereases when the density of base stations inereases. We willl leverage results from Q to 
ineorporate bloekage and differentiate between having LOS links and NLOS links from the base 
stations in the analysis of the problem of joint transmission in mmWave networks. 

In [Q (see also journal version in [1^) the authors used stoehastie geometry for studying 
mierowave joint transmission CoMP where single antenna base stations transmit the same data 
to single antenna users. Different performanee metries (ineluding eoverage probability) were 
eonsidered, to evaluate the performanee at the typieal user loeated at an arbitrary loeation (general 
user) and reeeiving data from base stations with the strongest average reeeived power. A user 
at the eell-eorner (worst-ease user) was also eonsidered. The eoverage probability was derived 
for both types of users under the assumption that the base stations have no CSI. The ease with 
full CSI was evaluated with different performanee metries (diversity gain and power gain). The 
derivation of the eoverage probability for a mmWave network with base station eooperation in this 
work is similar to that in Q for the general user, exeept that key faetors speeifie to the mmWave 
ehannel model have to ineorporated, some of whieh are the high direetional transmission at the 
base stations, bloekage and improved fading distribution due to sparse seattering. 

In Q the authors eonsidered the problem of finding a suitable single user MIMO transmit 
preeoding and reeeive eombining in mmWave systems under a set of hardware eonstraints 
suitable for large antenna arrays. Both problems (transmit preeoding and reeeiver eombining) 
were formulated as a sparsity eonstrained signal reeovery problem and solved using orthogonal 
matehing pursuit algorithms. The solution suggests that the transmitter applies a number of array 
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response veetors at the RF level (whieh are phase only veetors) and forms linear eombinations of 
these veetors using a digital preeoder. A similar observation was made for the reeeiver eombining 


operation. In [ 111 a multiuser MIMO downlink seheme, Joint Division Spatial Multiplexing 
(JSDM), was proposed. The seheme is suitable for frequeney division duplexing (FDD) systems 
with large number of antennas (massive MIMO) and non-ideal ehannel state information (CSI) 
at the base station. The base stations equipped with multiple antennas were assumed to serve 
K single-antenna users. Users who have identieal eovarianee matriees were grouped together 
while separate groups of users were assumed to have almost orthogonal eigenspaees of ehannel 
eovarianee matriees. The proposed two-stage JSDM exploits ehannel matriees properties and 
finds the optimal preeoding and prebeamforming matriees. In this paper, we assume that the 
eooperating base stations beam steer to the typieal user using veetors that ean be readily 
implemented using phase shifters, and the reeeiver applies a single veetor to proeess the reeeived 
signal from the eooperating base stations. 

Main Contributions: In this paper, we propose to study the benefits of base station 
eooperation in the downlink of a heterogenous mmWave eellular system as a mean to deerease 
signal outage. We antieipate that the eooperation provides substantial gain in eoverage with 
the antieipated improved fading distribution and extreme base station densifieation in mm Wave 
networks. Our extensive numerieal examples show that this is in faet the ease for the following 
seenarios. Case 1) for dense mmWave networks where the number of interfering LOS base 
stations inereases and Case 2) when there is no small seale fading ehannel on the ehannel gains 
from the eooperating base stations (a good assumption due to the high direetional transmission). 
We also provide examples when eooperation does not provide substantial inerease in eoverage 
probability. We eonsider a stoehastie geometry based model as in [|6|-p0|, p^ , to study eoverage 
in CoMP heterogenous mmWave network. To do so we need to ineorporate key faetors speeifie 
to a mmWave ehannel model. These speeifie mmWave eharaeteristies are: a realistie mmWave 
ehannel model, highly direetional ehannel gains and sensitivity to bloekages. 

Coverage probabilities are derived for the ease of a single antenna reeeiver (typieal user). We 
use eoneepts from [(^-[(^ to ineorporate bloekage, interferenee and different fading distributions 
(Rayleigh, Nagakami and no fading) in our analysis. The joint distribution of the eooperating 
base stations to the typieal user in the presenee of bloekage is also derived. 




5 


Paper Organization and Notations: The downlink CoMP mm Wave heterogenous network 
model, the beam steering at the base stations and the deeoding at the typical user are explained 
in Section The coverage probability in the absence of blockage, and with Rayleigh fading 
is derived in Section II-H In Section we consider Rayleigh fading mm Wave networks 
with a blockage parameter at each tier, and the coverage probability is derived accordingly. 
In Section we derive the coverage probability for the same network model with blockage 
but use the Nakagami fading distribution to model the fading distribution on the direct links of 
the cooperating base stations. The assumption of having no small scale fading for the channel 
gains from the cooperating base stations is further considered in Section |V| Proofs may be found 
in the Appendices. Tables |n IIT IV summarize all the notations used throughout the paper. 


TABLE I 

Poisson point process variables 


Notation 

Description 

K 

Total number of tiers 


Homogenous Poisson Point Process (PPP) indexed by k ^ [1 K] 

Afc 

Intensity of the PPP T>fc 

Pk 

Available power at each base station that belongs to tier k ^ [1 : K] 

V 

Points on 2D plane representing location of base stations 

IMI 

Distance from point v to the typical user located at the origin 

a 

Pathloss exponent assumed equal for all tiers 

Ofc = 6 4'ic} 

Normalized pathloss between each base station in and the typical user 

Afc(«) 

Intensity of 0 a; 

CD 

D 

II 

® 

Process representing the union of non-homogenous PPP, elements are indexed in increasing order WLOG 


Intensity of 0 

-Pj. 

Normalized pathloss 

II 

Set of normalized pathloss of the cooperating base stations 

/r'(7') 

Joint distribution of 7' 


II. Coverage Probability with no Blockage 

A. Network Model 

Consider a K tier heterogenous network where each tier is an independent two-dimensional 
homogenous Poisson point process (PPP). We denote the base station location process of tier 
k E [1 K] hy with density A^. The mmWave base stations that belong to the same tier 
k transmit with the same power Pk for k E [1 : K]. We study the coverage probability as 
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TABLE II 

General channel model variables 


Notation 

Description 

Nt,Nr 

Number of antennas at each base station and at the receiver 


MIMO channel from base station at location v to typical user 

hv 

Small scale fading 

L„ 

Number of channel clusters 

4>l 

Path angle at the transmitter 

4>l 

Path angle at the receiver 

fN 

Function that returns the index to which a base station at v belongs to 

_ 

- ||„||<. 

Pathloss 

^t(r) {■) 

Uniform linear airay vector representation at the transmitter (receiver) 

^t(r) 

Normalized transmit (receive) antenna separation 

Lt 

Normalized length of the transmit antenna array 

n 

Noise vector of i.i.d CA/^(0, crj) 


TABLE III 

Channel variables from the cooperating base stations 


Notation 

Description 

r 

Set of cooperating base stations with cardinality \T\ = n 

W, i6[l:|r|] 

Points on the 2D plane corresponding to cooperating base stations location (sometimes indexed by j instead of i) 


MIMO channel from the cooperating base stations 

hv- , 1 Ivi 

Channel parameters of the interfering links as defined in Table II 


Directional cosine given by 

x„, 

Transmit signal from cooperating base stations 


experienced by the typical user located at the origin, and denote the set of cooperating base 
stations, which jointly transmit to the typical user, by T G We assume that |T| = n, 

and that these n base-stations correspond to those with the strongest received power at the typical 
user receiver. In the rest of the section, we first describe the channel model and then derive the 
output signal at the typical user receiver. 


B. Simplified Clustered Channel Model 

A clustered channel model, [ [T3l , is used to model the wireless channel between the base 
stations and the typical user located at the origin. We assume all base stations have the same 
number of transmit antennas Nt, while the receiver has Nr receive antennas. The Nr x Nt channel 
matrix H.^, between a base station located at n G and the typical user is the sum of clusters 
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TABLE IV 

Channel variables from the interfering base stations 


Notation 

Description 

h, i G [1 : iri] 

Points on the 2D plane corresponding to interfering base stations locations 

Hi, 

MIMO channel from the interfering base stations 


Channel parameters of the interfering links as defined in Table II 


Directional cosine given by cos{(j)l.) 

Oi 

Angle used by interfering base station at position h to beam steer to a user other than typical user 

Qgt 

Directional cosine given by cos(0].) 

Xi, 

Transmit signal from interfering base stations 


and is expressed as 


where 


H, 




1=1 


( 1 ) 


p.. - 

• 'Iv,! = ii^ is the pathless, 

• f{v) is a funetion that returns the index k of the tier to whieh the base station at loeation 
V belongs to, 

. a is the pathless exponent, 

• IIuII is the distanee from the base station at loeation v to the user at the origin, 

• is the eomplex fading ehannel gain, 

• The veetors at(r)(0j)^^^) are the normalized uniform linear array (ULA) transmit and reeeive 
array response and are given by ^ Eq. (7.21), Eq. (7.25)] 

at(r)(0f)) = ^^[1, e-^'2^, • • • , (2) 

where A = 27rAt(r) cos(0t^^^) and At(r) is the normalized transmit (reeeive) antenna sepa¬ 
ration (normalized to the unit of the earner wavelength), at a path angle 0^ of departure 
(arrival) from the base station v. 

In the following, for simplieity, we shall eonsider the ease = 1. 
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A 



Fig. 1. A typical user is served by two cooperating base stations at locations v\ and V 2 , while being interfered by base station 
at location 


C. Received Signal at the Typical User 

In this section we will further divide the points n G into a set of points Vi and U to 
differentiate between the location of the cooperating and interfering base stations, respectively. 
The Nr X Nt desired channel matrices are denoted by H,,. for i G [1 : |T| = n], where n is a 
positive constant, while the interfering channel matrices are denoted by H;. for i G [1 : |T'^|]. 

Fig. shows an example of a network model, where two base stations at locations vi and V 2 , 
jointly transmit to the typical receiver located at the origin (indicated as Rxq) in the presence 
of a single interfering base station at location li. The MIMO channel matrices between the 
cooperating base stations and the typical user are given by and The channel matrix 
between the interfering base station and the typical user is denoted by The angles, 0*. 
and (pl., are the cluster’s angle of departure and arrival respectively from the base station Vi, 
i E [1 : 2], to the typical receiver. The angle (j)\^ is the angle of departure of the cluster from the 
interfering base station. The base station at li uses a beam steering angle 9i-^ to transmit data to 
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some other user (not the typical user) indicated as Rxi. The received signal is 

\T\=n in 

y= H.,X„. + y^H,.X,. + n 

i=l i=\ 

|T|.n \T‘\ 

= Y, + y] a,(4,) a.(4)*Xi, + n 


2=1 


2 = 1 


(3) 


where the first sum in Q is the desired signal from the mmWave cooperating base stations, 
while the second sum contains the signals from the interfering base stations. 

The user associates with a set of cooperating base stations T, that provide the strongest average 
received power as in p0| . Specifically, 

n p 

T = arg max (4) 

and := $fc\ T. The path angles, (/>*., i G [1 : |T|], and 0*,, i e [1 : \T^\], represent 
the angle of departure of the desired and interfering paths respectively, while i G [1 : |T|], 
and 0[,, i G [1 : |T^|], represent the angle of arrival of the received path from the cooperating 
and interfering base stations respectively. The transmit signals, X^., z G [1 ; |T|] and X;^, 
i G [1 : \T^\], represent the signal from the cooperating and interfering base stations within T 
and respectively, n is the noise vector of i.i.d C7\/^(0, components. 


D. Beam steering 

The base stations in T jointly send the same data to the receiver. Each base station beam 
steers to the typical user, therefore the transmitted signal is 

X„, = at(0*Js (5) 

for z G [1 : |T| = n], where s is channel input symbol transmitted by the cooperating base 
stations to the typical receiver. The signals transmitted by the interfering base stations are 

x,. = a.(»‘)s,, (6) 

for z G [1 : \T^\], where is the channel input symbol transmitted by the interfering base 
stations, while the angle 6\, is the angle used by base station Zj to beam steer to a user other than 
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the typical user, and is different from 0;. in general. We assume that s and si- are independent 
zero mean and unit variance random variables. 

Assumption 1: We assume that the cooperating base stations have perfectly beam steered to 
the typical receiver: notice that the angles in used by the base station to beam steer, are 
equal to the clusters’ angles of departure in the desired channel in 

E. Decoding 

The receiver uses a single vector w G to detect the scalar transmit symbol, that is, the 

processed received signal is given by 


V = w*y 

(7) 

n 

(8) 


i=i 


Remark 1. The choice of w in ([^ is one choice of a decoder that can be implemented readily 
using phase shifters in the RF domain (analog processing), in fact if one wants to consider a near 
optimal performance, then the work in Q, which finds a hybrid MIMO receiver combining algo¬ 
rithm and minimizes the mean-square-error between the transmitted and received signals under 
a set of RF hardware constraints for the resulting point-to-point channel should be generalized 
to finding a suitable algorithm for the downlink cooperative channel. 

Assumption 2: We assume perfect CSI of the path angles at the decoders since these angles vary 
slowly. However, we assume that the phases of the complex channel gains, G [1 ; |T|], 
are not available at the terminals as they change very quickly on the order of a wavelength and 
thus cannot be tracked. The performance here should be considered as an upper bound on the 
performance of the more realistic case with imperfect path angle assumption. 

F. Output Signal 

The output signal at the typical user under the previously stated assumptions is given by 

n / n 

y = w*y = ^ar(0^,)* I V NtNr 

j=i ^ i=i 

irq N 

i=l ' 


(9a) 
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= EE 

j=l i=l 
n |ri 

+ V V ^M.Gr{n^r, - % )si^ + ^ (%) 

f ^ ^ 3 ''i ‘i 

j=l i=l 

where 2; = w*n ~ CA^(0, a^), with cx^ = cr^w*w and where we introduced the antenna-array- 


gain functions 


G^iy) := sin{TTA N y) _ ^ 

iV^sin(7rA^I/) ' ^ 


( 10 ) 


^x(0i) (I02)) ^ ^ 

with f20 := cos(0) and A^;, x E {t,r} being the normalized antenna separation. 


( 11 ) 


G. SINK Expression 


Based on (|9bl), the instantaneous SINK is then given by 


SINK = 


X/ y/'yvihviCy 
i=l 


in 




i=l 


( 12 ) 


where G^. := Yl]=iGri^^r, - %i) and Di^ := 
Assuming a single antenna receiver with Nr 
in ([T^ simplifies to 


Yl]=l Gr{^4>l. 
- 1 (Q = Di^ 


n and = n^a‘^), the SINK 


n 

I y/'yvihvi I 

SINK =-’ (13) 

i + E -n,.W 

i=i 

The coverage probability for the typical user with SINK as in ( [T3| ) will be derived under 
the assumption that all angles are independent and uniformly distributed between [—ti.-Eti]. 
We will first assume that the receiver is present in a rich scattering environment (Rayleigh 
fading assumption), and in this scenario the coverage probability is given in Th. The case 
where each tier experiences blockage is then considered and the coverage probability is derived 
accordingly and is given in Th. The Nakagami fading distribution is then used to model 
the less severe fading distribution on the direct cooperating links and two upper bounds on the 
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coverage probability are then derived and are given in Th. and Corollary The ease where 
there is no small seale fading for the direet eooperating links is then eonsidered and the eoverage 
probability for this ease is given in Th. Future work ineludes deriving the eoverage probability 
for all the different eases deseribed above using ([T^, i.e, multiple antennas at the reeeivers. 


H. Performance Analysis 

Theorem 1. The coverage probability for the typical user, with a single antenna, in a downlink 
mmWave heterogenous network with K tiers, with base stations having ULA with Nt antennas, 
of which n jointly transmit to it, is given by 

P(SINR>T)= f (W) 

0<7^<---<7;,<+oo 

where 7 ' = ^ ^ : n], and the Laplace transform of the interference and the noise 




are given by 


i- r 

’1 n 

V A; 

J-2 \ 


]h{£) de 


X{v) dv ), 


(15) 

(16) 


C^(5)=e-*" 

where the antenna array gain Gt^e) is given by and the probability density function of 




/t(^) — 


fininll,!—£:} 


x{-i,-i-4 \7rV1 - {.£ + yY 
and fr'il') is the Joint distribution of y' = [ 7 ^, • • • , 7 ( 1 ] and is given by 


dy, 


(17) 


/r'(7') = 


-a(7;) 


(18) 


2 = 1 


while the intensity and intensity measure are given by 


K 


^—' a 

k=l 


(19) 


K 


A(7n) = '^T^XkPkl'd 


k=l 


Proof: Please refer to Appendix for the proof. 


( 20 ) 
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Fig. 2. Coverage probability in 03 for a two-tier network with parameters in Table |V] with two cooperating base stations 
(n = 2) and without base station cooperation (n = 1) and for different number of antennas. 


1. Numerical Results 

Example 1: In this section we numerically evaluate Th. We compute the coverage prob¬ 
ability for the typical user in a mmWave CoMP heterogenous network and compare it to the 
case with no base station cooperation. We consider a two tier network, K = 2 with parameters 
given in Table The noise variance is given by (T^(dBm) = —174 + 10 log^g(BW) + NF (dB), 
where BW and NF are abbreviations for bandwidth and noise figure respectively. In Fig. the 
coverage probability in ( [T?] ) for n = 2 and n = 1 is plotted. In the absence of blockage, the 
numerical results show that the increase in coverage probability with cooperation for the case of 
Nt = 8,16 antennas is almost 11% at T = 5 dB. While the increase is 10% for Nt = 32, 64 at 
T = 10 dB. As expected, an increase in the number of antennas at the base stations increases the 
coverage probability. For example, for the same threshold T = 10 dB, the coverage probability 
with cooperation and with Nt = 16 is approximately 0.5 while for Nt = 32 is 0.65. The increase 
in coverage probability can be interpreted as follows: the mmWave tier is relatively denser than 
that of a microwave tier, therefore the number of interfering base stations increases too. Thus, 
with cooperation limits interference and the coverage probability consequently increases. 

Remark 2. The authors in 0 compare a two tier network with parameters (power, noise. 
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TABLE V 


TABLE VI 


Tier for Eig. [^(Example 1) 


Tier parameters for Fig. [^(Example 2) 


Parameter 

Value 

Intensity 

Ai = (150V)"k A 2 = (50^)"^ 

Power 

Pi = 1 W (30 dBm) and P 2 = 0.25 W 

Path Loss 

a = 3 

Antennas 

Nt = 8,16, 32, 64 

Noise Figure (NF) 

10 dB 

Blockage 

Not Applicable 

Bandwidth 

1 GHz 


Parameter 

Value 

Intensity 

1 

(N 

0 

11 

Power 

Pi = 1 W 

Path Loss 

Ql = 2, £12 = 4 

Antennas 

Nt = 16 

Noise Figure 

5 dB 

Blockage 

P = 0.006,0.003,0.0143 

Bandwidth 

1 GHz 


intensities) suitable for microwave deployment with and without base station cooperation (n=2,3); 
an increase of 17 % was noted at a threshold T = 0 dB for the case of CoMP with two cooperating 
base stations when compared to the case of no cooperation. We shall show that a comparable 
gain (16%-18%) to the one reported in 0 can be attained with two cooperating base stations 
with non fading channel gains in Section |Vj Thus, the Rayleigh fading assumption considered 
here provides a worse case scenario. 


III. Coverage probability with Blockage 


A. Network model 


In this section we again consider a iC-tier heterogenous network where each tier is an 
independent two-dimensional homogenous (PPP). The base station location process of each tier 
is denoted by with density \k for k E [1 : K]. Each tier is characterized by a non-negative 
blockage constant {3^ for k G [1 : K] (determined by the density and average size of objects 
within the tier and where the average LOS range in a tier fc G [1 : iC] is consequently given 
by l//3k) as defined in [15| and used in Q, 0. Consequently, after defining the parameter 13k 
for k G [1 : K], we have that the probability of the communication link being a LOS link (no 
blockage on the link) within tier k is P(LOSfc) = 6“^'='', where r represents the length of the 
communication link, while the probability of a link being NLOS is P(NLOSa;) = 1 — P(LOSfc). 
The LOS and NLOS links will have different pathloss exponents, ai and 02 , respectively, and 
are the same for all k E [1 : K]. With the assumption of blockage the Laplace transform of the 
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interference in ( [T5] ) and the joint distribution of the cooperating base stations in ( [T8] ) have to be 
re-derived. 


B. Performance Analysis 

Theorem 2. The coverage probability for the typical user, with a single antenna, in a downlink 
mmWave heterogenous network with K tiers, and where each tier has a blockage parameter 
(3k, with n base stations having ULA with Nt antennas, jointly transmitting to it is given 
by ® but where now the intensity X{v) in ( [T5] ) is given by 


K 


X{v) = ^AkV‘^i ^ 


-afeu“l 




-bkv°‘2 


( 21 ) 


k=l 


_2_ _L _L _ 2 _ 

where Ak = nXk^PP , Uk = (3kPk ^, bk = (3kPk^, Bk = T^^k^Pk^, tind the distribution of the 


' <y.2 


distance of n closest base stations is given by (18) but where 


K 


A(7i) = E 


27rAi 






+ 7rAfc(7;Pfc)“2 


PI 


( 22 ) 


Proof: Please refer to Appendix for detailed proof. 


C. Numerical Results 

In this section we numerically evaluate Th. We compute the coverage probability for 
the typical user in a mm Wave CoMP heterogenous network in the presence of blockage. The 
examples provided illustrate scenarios when cooperation is beneficial (in terms of increasing 
the coverage probability) and examples when the increase is not substantial. Numerical results 
suggest that the former is in fact the case when the mmWave network is dense (captured 
by the tier radius and consequently its intensity) - a feature expected for millimeter wave 
networks 0. 0- This can be interpreted as follows, with extreme densification, the number 
of LOS interfering base stations increases and thus interference increases, a remark also noted 
in 0. Therefore, cooperation limits the interference by increasing the number of serving base 
stations and therefore providing higher coverage probabilities. 

Example 2: In Fig. we plot the coverage probabilities for a tier with parameters given 
in Table |V^ A tier with average radius of 80 meters is considered with blockage parameters 
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Fig. 3. Coverage probability in Th. |^for a one-tier network with parameters in Table |V^ with two cooperating base stations 
(n = 2) and without base station cooperation (n = 1) and for different blockage parameters. 


= 0.003,0.006,0.0143 (corresponding to average LOS range whieh is greater than 80 m 
for f3 = 0.003 and an average range that eannot reaeh a user at the edge for /9=0.0143). The 
coverage probability for this one tier CoMP mm Wave network with n = 2 and with base station 
power available P = IW is compared with the following cases, Case 1) a one tier mmWave 
network with no base station eooperation (n = 1), with a base station transmit power P = IW, 
Case 2) a one tier mmWave network with no base station eooperation n = 1, but with base 
station transmit power equal to the sum of transmit power if two base stations were to cooperate 
P = 2W. The inerease in eoverage probability for both oases is approximately an inorease of 0.12 
in probability for a threshold T = 5,10 dB. Moreover, it is interesting to note that an increase in 
the tier blookage parameter (shorter range of LOS links) would increase the coverage probability. 
This oan be interpreted as follows: an inorease in the blookage parameter inoreases the probability 
of blookage of the interfering LOS base stations, resulting in higher coverage probabilities. The 
curves oorresponding to the coverage probabilities to Case 1 and Case 2 are very olose sinoe 
the power at the interfering base stations has also inoreased with this assumption (whieh also 
means that this network is not noise limited). We shall show in the subsequent example, that this 
observation doesn’t hold for a less dense tier with a high probability of NLOS base stations. 
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Fig. 4. Coverage probability in Th. for a one-tier network with parameters in Table |VII| with two cooperating base stations 
(n = 2) and without base station cooperation (n = 1). 


Example 3: A tier with an average radius of 250 meters and with tier parameters given 


in Table VII and with a blockage parameter (3 = 0.02 (corresponding to a high probability 
of blockage and average LOS range of 50 m) is plotted in Fig. The increase in coverage 
probability due to cooperation in this case is minimal and is approximately 0.05 at all thresholds. 
This can be interpreted as follows: 1) a tier with high blockages will also block interfering signals 
and 2) when the density of base stations is not too dense, the n strongest base stations are not too 
strong to cause a substantial increase in coverage probability due to the fact that distance at which 
these cooperating base stations are located increases too (thus received power decreases). As seen 
in Fig. Q increasing the power at the base station (but no cooperation) provides higher coverage 
probability than the case with base station cooperation. We shall show that the observations made 
for this example do not hold when there is no fading on the direct links from the cooperating 
base stations in Section IVl 

Example 4: In an attempt to understand whether the observations hold for a network which 
is not noise limited as in Example 2 but for a larger number of antennas at the base stations, we 
consider the example of a dense mm Wave tier network (average radius of 50 m) and Nt = 64 
and with network parameters as in Table VIII In Fig. we plot the coverage probabilities 
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Fig. 5. Coverage probability in Th. |^for a one-tier network with parameters in Table [Vin| with two cooperating base stations 
(n = 2) and without base station cooperation (n = 1). 


TABLE VII 

Tier parameters for Fig. [^(Example 3) 


Parameter 

Value 

Intensity 

A = (250V)“^ 

Power 

Pi = 1 W 

Path Loss 

ai = 2, a2 = 4 

Antennas 

II 

Noise Figure 

5 dB 

Blockage 

/3 = 0.02 

Bandwidth 

1 GHz 


TABLE VIII 

Tier parameters for Fig. [^(Example 4) 


Parameter 

Value 

Intensity 

II 

o 

to 

1 

Power 

Pi = 1 w 

Path Loss 

II 

Ol 

Q 

cL 

II 

s 

Antennas 

Nt = 64 

Noise Figure 

5 dB 

Blockage 

P = 0.004 

Bandwidth 

1 GHz 


corresponding for the different eases whieh are deseribed in Example 1. The observations made 
in Example 1 hold for this example too with almost the same inerease (11%) in eoverage 
probability for a threshold T = 10,15 dB. 

IV. Coverage Probability with Nakagami fading and blockage 

In this seetion we eonsider the same network model as in Section IlII-AI hut choose a different 
fading distribution on the ehannel gains from the strongest eooperating base stations, similar 
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to Q. In particular we consider Nakagami fading with parameter m, while keeping the same 
assumption of Rayleigh fading for the interfering channel gains. 


Using the coverage probability expression for a general fading distribution [16, Eq. 2.11], 
we are then able to derive an upper bound on the coverage probability for this network. We 
then consider another upper bound by evaluating the network in the absence of interference. 
The coverage probability is then defined as the probability that the signal-to-noise-ratio (SNR) 
is greater than a certain threshold. Using complex analysis methods of integration we derive 
closed form for the coverage probability. 


A. Performance Analysis 

Theorem 3. An upper bound on the coverage probability for the typical user, in a downlink 
mmWave heterogenous network with blockage with K tiers, where each tier has a blockage 
parameter /3k, and with n base stations having ULA with Nt antennas jointly transmitting to it, 
with the assumption of Nakagami fading on the cooperating channel gains is 


P(SINR > T) < 


friy) [ Ci{2j7;T's)CN{‘2j7iT' ^ - -ds dy' 


where 7 ' = 


17^11 


0<7]<---<7^<+oo 


for i G [1 : n] and T' = 


2jTTS 


(23) 


7i^ 


while the joint distribution of 7' = 


bi) ■ ■ ■ given by ( [T8] ) and where the Laplace transform of I (assuming Rayleigh fading 

on the interfering links) is given by while the intensity and intensity measure are given 


by (21) and (22) respectively, and where 


/-UP 


's) = 


(1 + s/m)^'^ 

C^(s) = 

Proof: Please refer to Appendix for a detailed proof. 


(24) 

(25) 


Remark 3. For the case of no base station cooperation n = 1, the coverage probability in ( [23] ) 
above is exact and is not an upper bound. 

Remark 4. If one desires the exact coverage probability, one can obtain it with the theorem 

n 

below, with replaced by the true Laplace transform of the signal S = | ^ 


i=l 
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TABLE IX TABLE X 

Tier Parameters for Eig[^(Example 5) Tier parameters for Eig.|7(a)|(Example 6) 


Parameter 

Value 

Intensity 

Ai = (200V)“^ 

Power 

Pi = 1 W (30 dBm) 

Path Loss 

ai = 2 and 02 = 4 

Antennas 

Nt = 16 

Noise Figure 

5 dB 

Bandwidth 

1 GHz 

Blockage 

0.025 


Parameter 

Value 

Intensity 

> 

0 

0 

to 

1 

Power 

Pi = 1 w 

Path Loss 

ai = 2, 02 = 4 

Antennas 

11 

Noise Figure 

5 dB 

Blockage 

P = 0.025 

Bandwidth 

1 GHz 


Corollary 4. An upper bound on the coverage probability in the absence of interference for the 
typical user, in a downlink mmWave heterogenous network with blockage with K tiers, where 
each tier has a blockage parameter jdk, and with n base stations having ULA with Nt antennas 
jointly transmitting to it, with the assumption of Nakagami fading on the cooperating channel 
gains is 

P(SNR > T) < J friy') dy'. (26) 

0<7j<---<7;<+oo 


where z* = a- = 

vviicic 27 rj, h p 


for i G [1 : n] and T' = ^ ,_i while the joint distribution of 


y = b'l,--- ,yn] is given by G1 with an intensity and intensity measure as in ( [2T| ) and 
respectively and where 


9b) = {-If 


1 _ n — , 2 


(^)nm(2vrj2;) 

and where is the {nm — 1) derivative of the function g{z). 


(27) 


Proof: Please refer to Appendix for the proof. ■ 

A similar remark to that made in Remark 4 ean be made for the eoverage probability (in 


absenee of interferenee) in (26). 


B. Numerical Results 

Example 5: We eonsider a one tier network with two eooperating base stations n = 2 and 


with tier parameters given in Table IX The eoverage probability in ([2^ for the ease of m = 3 
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Fig. 6. Upperbounds on coverage probability in Th. |^and Corollary for a one tier network with parameters in Table [I^ and 
for n = 2 (two cooperating base stations) and for n = 1 (no cooperation). The lower bounds are plotted using Th. 


with and without base station cooperation are plotted. The purpose of this numerical example 
is to show that for tiers with high probability of blockage, in this case taken to be /3 = 0.025 
(corresponding to a high probability of blockage and average LOS range of 40 m), evaluation of 
the coverage probability of the network with and with the absence of interference yields almost 
exact numerical results. Therefore, we fix n = 2 and we plot the coverage probability in ( |2^ 
and ( [^ for Case 1) m = 3 and n = 2. While we use Th. to plot the coverage probability 
for the Rayleigh fading Case 2) m = 1 and n = 2 (also with and without interference). The 
two curves shown in Fig. |^for each of the cases corresponding to the Rayleigh and Nakagami 
fading almost exactly overlap. 

V. Coverage Probability with No Small Scale Fading 

In this section we consider the same network model as in Section ITTT-AI but where the 
cooperating channel gains do not experience any fading. As shown in 0, the assumption of 
having no small scale fading from the serving base stations is a good assumption in mmWave 
systems due to the highly directional transmission and when the receivers are not present in a rich 
scattering environment (in rich scattering environments Rayleigh fading may be more reasonable). 
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The Rayleigh distribution is used to model the fading distribution of the interfering ehannel gains. 
Interestingly, we will show through numerieal examples that the inerease in eoverage probability 
with two eooperating base stations is more pronounced when the cooperating channels experience 
no fading than that obtained when the fading is assumed to be the Rayleigh fading. 


Theorem 5. The coverage probability in the absence of small scale fading for the typical user, 
in a downlink mmWave heterogenous network with blockage with K tiers, where each tier 
has a blockage parameter 13k, and with n base stations having ULA with Nt antennas jointly 
transmitting to it, with the assumption of no fading on the cooperating channel gains is 

P(SINR>T)= [ frh') [ Ci{2jTiTs)CN{2jnTs) ^^^~l^''"^ ~ ^ ds dy' 

J J—CO ZjTTS 

0<7j<-"<7^^<+oo 


(28) 


where 7 ' = for i G [1 : n] while the joint distribution of y' = [ 7 ^, • • • is given 

by ( [T8] ) with an intensity and intensity measure as in pT| ) and ( [22] ) respectively and where 
l^s{s) = and 


Proof: The proof follows easily from finding the Laplace transform of the desired signal. 


A. Numerical Results 

Example 6 : In this example we numerically evaluate Th. [^ We consider a one tier network 
with tier parameters given in Table [^ The coverage probability in Th. [^ with two cooperating 
base stations n = 2 is plotted against the two different cases described as in Example 2 of 


Section III-B We also plot the three curves corresponding to the Rayleigh fading assumption 
on the cooperating channel gains in Th The curves in Fig. 7(a) suggest that the increase in 
coverage probability with CoMP is mainly due to a power increase since the curve corresponding 
to the case when there is no cooperation but with double transmit power is very close to that 
of the curve of having two cooperating base stations but with half the transmit power. This also 
implies that this network is noise limited. This observation is valid for the networks where gains 
from cooperating base stations are non fading or when they are Rayleigh fading. It is interesting 
to note that the increase in probability with two cooperating base stations (almost 18 % for 
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(a) Coverage probability in Th. for a one tier network 
with parameters in Table with two cooperating base 
stations (n = 2) and with no cooperation (n = 1). The 
coverage probabilities for the same network parameters 
but with Rayleigh assumption in Th. are also plotted 
for (n = 1) and (n = 2). 



(b) Coverage probability in Th. for a one tier network 
with parameters in Table |VII| for j3 — 0.003 with two 
cooperating base stations (n = 2) and with no cooperation 
(n = 1). The coverage probabilities for the same network 
parameters but with Rayleigh assumption in Th.[^are also 
plotted for (n = 1) and (n = 2). 


threshold T = 5 dB) for the ease when there is no small seale fading exeeeds the inerease with 
CoMP (5 % for the same threshold) for the ease when the channel gains are Rayleigh distributed 
and hence implying that CoMP provides larger gains for mmWave networks with no small scale 
fading from the serving base stations. 

Example 7: In an attempt to understand whether the observations made in Example 6 hold 
for a network which is not noise limited, we consider the network in Example 2 (/3=0.003) and 


compare the Rayleigh fading case to the case when there is no small scale fading. In Eig. 7(b) 


we plot the curves corresponding to the different cases as explained in Example 2 along with the 
curves corresponding to the coverage probabilities non fading networks in Th. The increase 
in coverage probability with CoMP {n = 2) is an increase of 16% at threshold T = 10, while 
the increase with CoMP for the Rayleigh fading case is almost 12 % at the same threshold. 
Therefore, the observation made in Example 6 does in fact hold for this network too. 


VI. Conclusions 

In this paper we have considered the problem of base station cooperation in mmWave het¬ 
erogenous networks. Using stochastic geometry, coverage probabilities were derived at the typical 
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user, accounting for directionality at the base stations, blockage, interference and different fading 
distributions (Rayleigh and Nakagami). Numerical results suggest that coverage with cooperation 
rival that with no cooperation especially in dense mmWave networks with no small scale fading 
on the cooperating channel gains. Future work includes deriving the coverage probability at a 
multi-antenna typical receiver and accounting for possible errors due to beam steering. 


Appendix A 
Proof of Theorem [H 

The analysis of the coverage probability for the mmWave heterogeneous network is similar 
to that in [|^ Appendix A] with two major differences. The first difference is the presence of 
multiple antennas at the transmitter. The second difference is that the interference is a function 
of i.i.d uniformly distributed random variables, assuming that the path angles are independent 
and uniformly distributed over [—tt, +7r]. 

Let 0fc = G for fc e [1 : iF]. Its density can be derived using the Mapping 

theorem Thm. 2.34] and is given by 

OtT 1 2 

Afc(n) = Xk—PkV~^, fc G [1 : K]. (29) 

The process 0 = U^^0fc has a density 

K 

\{v) = Y,^k{v). (30) 

k=l 

A. Distribution of Strongest Base stations 


We assume that the elements in the process 0 are indexed in increasing order. Let 


li = 


P 


/K) 


then 7 ' = {yj, • ■ ■ , 7 ^^} denotes the set of normalized pathloss of the cooperating base stations. 
We first present the distribution of the two nearest base stations by following similar steps as 


done in [ [r 8 | , then derive the distribution of n closest base stations. The distribution of the closest 
two base stations (assuming two cooperating base stations) is given by 


/r'(71,72) - /r'|ri( 72 l 7 i)/r'( 7 i) 


( 31 ) 
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where the distribution of the first base station with strong reeeived power is obtained from the 
null probability of a PPP is 

/ri(7',)=A(7:)e-''«' (32) 


while the eonditional distribution is given by 

/rilr;(7il7:) =A(7i)e-''<’;>+^«' 


(33) 


The joint distribution for the ease of n = 2 base stations is obtained by substituting ( |32| ) and ( |33] ) 
in ( |3T] ). The result ean be generalized to any number n of eooperating base stations 


fr'W) 


-a(7;) 


(34) 


2=1 


B. Derivation of Coverage Probability 

We have assumed that the eooperating base stations have normalized pathloss ii by i < n, 
then the desired signal power at the numerator of ( [T3| ) ean be re-written as 

n 

2=1 i<n 

We have that the interfering base stations are indexed with i > n, then the power of the 
interferenee / ean be expressed (by replaeing the index /j with just i) as 

irr 


I = 7jfti.nG,((i«| - Si,;)!- 

• ^ 4 4 

2=1 

= Z2r‘|fc.nG,(Ti)p 


(35) 


i>n 


The eoverage probability for a threshold T, ean be re-written as 


a 


P(SINR > T) = P ( 5 > T(/ + —) ) = Eyj 


i*(lE 34 %r>r(/ + ^ 


i<n 


= Ey,/ 


(c) 


exp 


'^i<n T 


2 . 
Nf 

/-I 


= 




C, 


T 




n \^N 


T 


C 


T 


N 


'Thikn T 


/-I 


'^i<n T 


73T MV 


0<7(<"-<7^j<+oo 

where (a) follows from the eumulative density funetion of the exponentially distributed random 
variable S (due to Rayleigh fading assumption) with mean Yl (b) follows from the definition 

i<n 
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of the Laplace transform of I, Ci{s) = E[e and the Laplace transform of the noise, Cn{s) = 
(c) by definition of the expectation with respect to the distribution of 7'. 

Next we evaluate the Laplace transform of the interference I, but before going into the details 
of the derivation we need to find the distribution of T, := , since the interference 

in is a function of the beam forming gain function which in turn is a function of Tj. The 
beam forming gain is given by dTO] ). 

With the assumption that the interfering path angles of departure and the beam steering angle 
used by the interfering base stations are i.i.d ~ [/([—tt,+ 7 r]), then the directional cosine Vt^t 
and Qgt are random variables with the following common probability density function 


fn{uj) = { 


if —1 < a; < 1; 


0 otherwise. 

then the distribution of Tj = — VLqi is the result of the convolution of the probability 


density functions of and fig* and is given by 

^ 7 . 


h,{£i) - 


r’minjl,!— 


,x{-i-i-£7 - {et + \/l - 

Then the Laplace transform of the interference can be derived 


dy 


Ci{s) = E 

g i>n 

= E 





_ i>n ^ _ 


E 


{Ti},e 


= E{Tj,e 


He 


\h\ 


t>n 


2 e 




n 


i>n 


= Ee 


= E0 


He, 

. i>n 

n 


1 + o|G,(T,)| 27 '-‘ 
1 

r-+2 / 


i>n 


-2 \ 1 +'S|Gt(£)p 7 ; 


731 ]h{£) de 


(e) 

= exp I - 


hn 


r+2 


1 - 


^ / t (£) de 


X{v) dv 


( 36 ) 


( 37 ) 


1-2 y 1 +'S|Gt(£) pu" 

where (a) follows from the i.i.d distribution of and their independence from 0 and Tg 
where (b) follows from the Rayleigh fading assumption and the moment generating function of 
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an exponential random variable; where (e) follows from the i.i.d distribution of Tj and their 
independenee from 0; (d) from the taking the expeetation with respeet to the random variable 


Tj whose distribution is given by ( |3^ ; (e) follows from the probability generating funetion of 
poisson point proeess |[^ Thm. 4.9] (as used in [ 10 Eq. (38)]) and where X{v) is given by pO] ). 

Next we give an approximation of the Laplaee transform of the interferenee for easier nu- 
merieal evaluations (by approximating the beam forming gain funetion by a pieeewise linear 


funetion) and to eompare our results with p0[ , the Laplaee transform of the interferenee is 

-+2 


Cl is) Ee 



I + s\Gtie)\^ii 


731 /t(£) de 


-l/Lt fl/L 

frie) de + 

2 


I — - j^fr(e)de+ hie) de 

-l/Lt 1 + S7i Jl/Lt 


c S 7 , 


-1 


-1 


. i>n 


ia) 

= exp - 


c sv 


-1 


1 + su 


-1 


A(u) dv 


1 + sYi 

where in (e) an approximation of the gain funetion was used whieh is given 

I 1 -7- < E < f , L, = N,A, 

G,(£)= ] 

I 0 if otherwise. 

(f) defining c := hie) de. 

Remark 5. If c = 1 then the Laplaee transform in step (g) simplifies to that in p0[ Eq. (38)]. 


Appendix B 
Proof of Theorem 


A. Intensity and Intensity Measure 


Let 0fc = V e for k E [1 : K] with intensity Xkiv) given in ( |29l ). The pathloss a is a 


random variable that takes on values ai and 02 with probability e and 1 




respeetively 


(note that we have dropped the {|.|| of v for easier notation). Then the proeess 0 = Uhhk 
is a non-homogenous PPP with density A(u) = Ylk=i hiv). In the following we eompute the 
intensity and intensity measure of 0^ for k E [1 : iT]. By using the Mapping Theorem p7| 
Thm. 2.34] the intensity measure and the intensity of eaeh tier fc, fc G [1 : A'], are given by 






2TrXkVe + 


qrPfe)“2 


27rAfcr;(l — e ^'‘^)dv 
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2nXk 

PI 

27r 

~w 


1-e (1+ /gfc(rPfc)“i)^ +7rAfe(rPfc)“2 

^l_e-/^'=(^^'=)^(l + /?fc(rPfc)^) 


(iu 

1 


with Ak = TrXk^Pk^ > , bk = f^kPP'^ and Bk = TrXk^PP ''. 


' Q 2 


The process 0 = U^^^Qk has the following intensity measure and intensity 


K 


Hv) = 

k=l 

K K 

H^) = ^kiv) = ^ AkV 


P-1 ^-a^.v°‘l 


“1 e 




— b^v “2 


k=l 


k=l 


(38) 

(39) 


(40) 

(41) 


Appendix C 

Proof of Theorem [3] and Corollary [4] 

A. Proof of Theorem 

Let us re-consider a different distribution on the direct links - while keeping the same Rayleigh 
fading assumption on the interfering links - in particular let us consider that the fading is 
Nakagami with shape parameter m and scale parameter 6^ = 1. In this case we will derive 
the distribution of an upper bound on the desired signal in particular the distribution of 

_ ^UP 

i<n i<n i<n 

Then we have 

P (s“'’ > T{I + ^)) = p [e + W,' 


s=iEP'P|PE2r'Ei'".i^ = s' 


, i<n 


\T'=-P^,svp=j:,\h 

Ei 7i 


The distribution of the upper bound and its Laplace transform are respectively given by 

1 


El'* i I ^ Gamma(nm, 1/m) -fA £ 5 up(s) 


i<n 


(1 + s/m)' 


We then have from [16, Eq. 2.11] that the coverage probability is 

( 2 
S > 
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= [ fr'i'y') [ Ci{2j7TT's)CN{2j7rT's) ^^''^^ - -ds di 

J J —oo ZJTTS 

0<7j<-"<7^<+oo 

where the joint distribution of 7 ' is given by 

n 

i=l 

Next we have from ( [TT] ) the Laplaee transform of the interferenee Ci{s) with an intensity X{v) 

given by ( |4T] ) while the Laplaee transform of the noise is given by 

2 2 

Cn{s) = E[e"^^] = (43) 


B. Proof of Corollary 


The eoverage probability in the absenee of interferenee is given by ( [23] ) with Ci{s) = 1 is 


P(SNR>T)= j frfi) 

0 < 7 (<---< 7 ^<+oo 


' m. ' 


- 1 


X) 2j7rs 

poo 

/r'( 7 ') / f{s)dsdy'. 


3 Nt ds dy' 


(44a) 

(44b) 


0<7(<---<7;,<+oo 


-V-" 

Q 


In the following we seek to solve Q. Note that a pole of order nm exists in the integrand thus 
the integral Q ean be solved using eontour integration and its Residue (Res) is 


nm—1 


Q = Res,*-^ [f(z)] = lim --^^ ) (z- z*r^f(z). 

V {nm-l)\\dzj ^ ! 

The funetion f{z) ean be re-written in the following form 


(45) 


f{z) = 


9{z) 


{-If 


1 /-I 27 rjz \nm 

'' ^ _ m ' 

(27rj)’*"* {2-Kjz) 


e 


(z — z*)’’ 


9 {z) = {-ly 


7 - (1 - 27rjz)-- 


(27rj)"-™-(27rjz) 

Then after substituting the funetions in ( |45| ), we ean express the integral I as 


(46a) 

(46b) 


Q = 


{nm — 1 )! 


(47) 
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